Ground-state properties of hard-core anyons in one-dimensional optical lattices 
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We investigate the ground-state properties of anyons confined in one-dimensional optical lattices 
with a weak harmonic trap using the exact numerical method based on Jordan- Wigner transforma- 
tion. It is shown that in the Bose limit (x = 1) and Fermi limit (x = 0) the momentum distributions 
are symmetric but in between they are asymmetric. It turns out that the origin of asymmetry comes 
from the fractional statistics that anyons obey. The occupation distribution and the modulus of 
natural orbitals show crossover behaviors from the Bose limit to the Fermi limit. 
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I. INTRODUCTION 

The physical systems with fractional statistics have 
been a subject of great interest in past decades and 
been intensively studied for two-dimensional systems 
[H, 0, H, 0] ■ For instance, the elementary excitations of a 
fractional quantum Hall (FQH) liquid satisfy fractional 
statistics. This has been observed in two-dimensional 
electron gas and anyon has become an important con- 
cept in studying the FQH effect [H, H, @ . As a natural 
generalization of the Bose and Fermi gas, anyon gas has 
also found application in various one-dimensional (ID) 
systems [1, 0, H, [t|. Despite no explicitly experimen- 
tal proof of the realization of ID anyon gas, many the- 
orists have dedicated to study them @, H, El EL EH 
El, El, El, El, El El El M H3. Currently anyons 
also stimulated intensive research on topological quan- 
tum computation because the statistical properties are 
closely related to the topological order. Particularly, by 
controlling exchange interaction between pairs of neutral 
atoms in optical lattice, Anderlini et al. [22| realized the 
key operation in quantum information processing, i.e., 
the quantum SWAP gate (which control the states in- 
terchange between two qubits). Cold atoms in optical 
lattice are also proposed to create, manipulate, and test 
anyons [H, H [25[ . 

By tightly confining the particle motion in two direc- 
tions to zero point oscillations ID quantum gas is ob- 
tained [H, [27], Hl| , where the radial degrees of freedom 
are frozen and the quantum gas is effectively described by 
a ID model [29]. Experimentally, by means of anisotropic 
magnetic trap or two-dimensional optical lattice poten- 
tials, a ID Bose gas in the strongly correlated Tonks- 
Girardeau (TG) regime can be achieved [HEH. The TG 
gas has been shown to display the " fermionized" char- 
acter in many aspects. For example, it has the same 
density distribution and thermodynamic behavior as the 
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free fermion [iHU, IH, IHIl^l . However, the off-diagonal 
density matrices and the momentum distributions exhibit 
quite different behaviors due to the different permuta- 
tion symmetries of the Bose and Fermi wave functions 
[351 l3a ] . The momentum distribution of fermions shows 
typical oscillations but that of bosons shows the struc- 
ture of single peak. Both of them are symmetric about 
the zero momentum. For the anyon gas satisfying frac- 
tional statistics, however the momentum distribution is 
shown to be asymmetric when the statistical parameters 
deviate from the Bose and Fermi limit [H, El, El E^ , 
and with the change in statistical parameter the distri- 
bution evolves from a Bose distribution to a Fermi one 
and vice versa. 

While most of the theoretical works on the ID anyon 
gas focus on the continuum system, in this paper we in- 
vestigate the ground state of a ID anyonic system con- 
fined in optical lattice with weak harmonic trap in the 
hard core limit. Although the hard-core bosons have been 
studied extensively and intensively for both homogeneous 
continuum system and lattice system, no result was given 
for the hard-core anyons (HCAs) in optical lattice. We 
shall study how the fractional statistics affect the ground- 
state properties, such as the momentum distribution. By 
extending the exact numerical method originally used to 
treat hard-core bosons by Rigol and Muramatsu [3?J to 
deal with the hard-core anyons, we evaluate the exact 
one-particle Green's function of the ground state, with 
which the reduced one body density matrix (ROBDM) 
and thus the momentum distribution can be obtained 
exactly for different statistical parameters. This method 
is based on Jordan- Wigner transformation and has been 
applied in the investigation of hard core bosons confined 
in optical lattice for both ground state and dynamics. 
It turns out to be very efficient to study the universal 
behaviors of the system with arbitrary confining poten- 
tials combined with one-dimensional optical lattices. The 
properties of anyonic statistics shall be displayed in the 
momentum distribution and with the change in statistics 
parameter the system exhibits the Bose statistics, Fermi 
statistics, and the fractional statistics in between. The 
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mathematical origin of the asymmetry can be traced back 
to the reduced one body density matrix. 

The paper is organized as follows. In Sec. II, we give 
a brief review of ID anyonic model and introduce the 
numerical method. In Sec. Ill, we present the momen- 
tum distributions, ROBDM, and occupation for different 
statistics parameter and filling numbers. A brief sum- 
mary is given in Sec. IV. 



II. FORMULATION OF THE MODEL AND 
METHOD 

The second quantized Hamiltonian of the hard-core 
anyons confined in optical lattice of L sites with a weak 
harmonic trap takes the form of 
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with the lattice dependent external potential 

Vi = Vb(/-(L + l)/2) 2 , 

where Vq denotes the strength of the harmonic trap. The 
anyonic creation operator a\ and annihilation operator Oj 
satisfy the generalized commutation relations 



a, a 



-iX-xe(j-l) f 
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a jai = -e lx ^-^a iaj (2) 

for j 7^ I with the addition of hard-core condition af — 

a\ 2 = and |a;,aj| = 1. The sign function e(x) gives 

-1, 0, or 1 depending on whether x is negative, zero, 
or positive. The parameter \ is related with fractional 
statistics and will be restricted in the regime of [0, 1] in 
the present paper. Particularly, \ — an( l 1 correspond 
to Fermi statistics and Bose statistics, respectively. In 
the Hamiltonian the particle number operator m — a) at, 
and t denotes the hopping between the nearest neighbor 
sites. 

This model can be solved exactly by the generalized 
Jordan- Wigner transformation, 

a 3 = exp i X n f»f» fi> 
V !< s <i / 



a] = fj exp j -ixTr 
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where /• and fj are creation and annihilation operators 
for spinless fermions. The hard-core anyonic Hamiltonian 
with N anyons can be mapped onto the noninteracting 
fermionic system for Np fermions (Np = N), 



H F = -tj2 (fi+Ji + H - c -) + E v ™ w 



i=i 



with fermionic particle number operator m — ////.. 
While the eigen problem for Hamiltonian Hp for van- 
ishing harmonic potential Vi = can be obtained easily 
through Fourier transformation, we investigate here the 
situation with weak harmonic trap and Rigol-Muramatsu 
method [37| should be a good choice. Using the above 
transformation the one-particle Green's function of hard- 
core anyon is formulated as 
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with 
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I^hca) 1S the ground state of hard core anyonic system 
and |*^) is the ground state of free spinless fermionic 
system. The Green's function can be obtained by con- 
structing the many-particle ground state of fermions with 
the eigenstates of single-particle, 

|a> = 4|0>=5>° (0//|0), 
i 

where a means the ath state and I means the Zth site. 
The many body ground state of Nf free spinless Fermions 
takes the following form: 
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with Pi n = ip n (I), which can be expressed as an L x Nf 
matrix P. After an easy evaluation the state \^p) reads 



Nf+l l 



71=1 1 = 1 



with 



Pin = exp (-iXTT) Pl„ 
Pin 
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r ln 



for I < j — 1 
for I > j 



for n < Nt, and P 



lA 



i f , ana f^y+i = 1 and P$ f+1 =0(1^ j)- The 
state \^p) has the same form with the replacement of j 
by I. The Green's function is a determinant dependent 
on the L x (Nf + 1) matrices P /A and P' B , 
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FIG. 1: (Color online) Momentum distribution of the ground 
state for Vo = 0, N = 50, and L = 300. The unit of k is l/o 
with a the lattice constant. For convenience, hereafter we set 

0=1. 



Using the anyonic commutation relation the ROBDM can 
be expressed as 



/ = (a] ai ) = 5 fl (1 - G jt ) - (1 - 6 jl )e- i x*G jl . (8) 



Pi 



The natural orbitals <$P are defined as the eigenfunctions 
of the one-particle density matrix, 



3=1 



(9) 



and can be understood as the effective single-particle 
states with occupations A,,. In order to investigate the 
effect of statistical parameter \ we Wlu focus on the mo- 
mentum distribution, which is defined as 



(10) 
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III. PROPERTIES OF THE GROUND STATE 

We first investigate the momentum distributions for 
the situation in the absence of a harmonic trap. In this 
case the single-particle eigenstate for t = 1.0 can be for- 
mulated as (f>ctiP) — sin cxItt / -\- 1 ) / ■s/ C a with the nor- 
malized constant C a = ^,- =1 sin 2 ajw/(L + 1) for the 
fixed boundary condition [38| . Following the procedure 
in Sec. II the momentum distribution can be obtained 
exactly and is shown in Fig. 1. The similar behavior 
as the finite continuum model is manifested (l6j . In the 
Bose limit (x = 1) most of bosons populate at the zero 
momentum state, which corresponds to a Bose super- 
fluid, while in the Fermi limit (x = 0) the step-function 



FIG. 2: (Color online) Momentum distribution of the ground 
state for V = 1.0 x 10~H, N = 50 (a) and (b) 150, and 
L = 300. 



distribution is shown, which is the characteristic feature 
of free fermions. Both of them are symmetric about the 
zero momentum. When the statistical parameter \ de- 
viates from these two limits the momentum distribution 
is asymmetric about the zero momentum. Anyons dis- 
tribute in the regime of positive momentum with more 
probability than that in the regime of negative momen- 
tum and redistribute between these two regimes with the 
increase in statistical parameter \. Finally the distribu- 
tion evolves between Fermi and Bose distribution with 
the change in statistical parameter. 

When the harmonic trap is present we have to turn to 
the numerical diagonal method [13] in order to get the 
N lowest single particle eigenstates of free spinless Fermi 
model and construct the ground state of many body any- 
onic system. The momentum distributions are displayed 
in Fig. 2 for filling number N — 50 and N = 150 confined 
in an optical lattice with L = 300 sites. In this situation 
the anyons populate in broader regime compared with the 
confinement- free case as shown in Fig. 1. This is more 
evidently seen in the Fermi limit. In the absence of har- 
monic trap anyons (pure Fermions in this limit) occupy 
the momentum states in the regime of— 1/2 < k < 1/2 
almost homogeneously while the combined harmonic trap 
results that higher momentum states are occupied and 
the distribution is no longer a sharp step-function. For 
all statistical parameters the half width of momentum 
distribution becomes wider. This can be understood by 
uncertainty principle that the presence of harmonic trap 
reduces the uncertainty in the coordinate space and leads 
to the increase of uncertainty in the momentum space. 
When more particles are loaded in the lattice they will 
populate at higher momentum states with more proba- 
bility and the peak of momentum profiles will not be so 
sharp as those in the case of less filling number. Ob- 
viously the filling number and the strength of harmonic 
trap shall not affect the statistical property that with the 
change in statical parameter anyon's momentum distri- 
bution evolves continuously from a Bose distribution to 
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a Fermi one. In these two limits the distributions are 
symmetry about the zero momentum and in between the 
profiles are asymmetric. 

In order to clarify the origin of asymmetric momentum 
distribution of anyons we show the ROBDM for L = 300 
lattice sites with 50 anyons in Fig. 3, where j and I 
ranged from 1 to L. In the Bose and Fermi limits the 
ROBDMs are real (top row) whereas in the regime devi- 
ating from these two limits the ROBDMs are Hermitian, 
i.e., possessing symmetric real part (Re[pji]=Re[pij]) and 
anti-symmetric imaginary part (Im[/Oj2]=-Im[/9jj]), which 
are exhibited in last four rows for \ =0.2, 0.4, 0.6, and 
0.8. Therefore the formula of momentum distribution 
can be departed into two parts as below that one part is 
even function of k and the other one is odd function of 
k, 



77 3,1=1 
1 L 

= — ^ {R- e [Pjl] cos k (j - I) + Im [p i] sin k(j -I)}. 
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It is because of this that the momentum distribution be- 
comes asymmetric about zero momentum. The imagi- 
nary part Impji is an odd function of statistical parame- 
ter x so the peak at positive momentum as shown above 
will shift to negative momentum if we take x as negative 

(x--x)Q. 

In Fig. 4a we display the occupation of natural orbitals 
for the system with 50 anyons in 300 lattice sites com- 
bined with weak harmonic trap (Vq = l.Ox 10~ 4 £). In the 
Fermi limit each anyon occupies one orbital and the oc- 
cupation distribution is a step function. In the Bose limit 
most anyons occupy the lower orbitals and the occupa- 
tion distribution exhibits sharp single-peak structure. In 
the regime deviating from these two extreme points the 
evolution from one to the other indicates that with the 
decrease in statistical parameter more and more anyons 
occupy higher orbitals and finally distribute in the low- 
est N orbitals homogeneously in the Fermi limit. The 
statistical effect on the natural orbital is shown in Figs. 
4b-4d for the same system as Fig. 4a, where the modulus 
of the lowest natural orbitals ({(j) 1 ]) is exhibited for sta- 
tistical parameter x=0.0, 0.4, and 1.0. In the Bose limit 
the modulus are almost smooth and only weak oscilla- 
tions appear around the boundary while the oscillations 
become more and more obvious with the decrease in sta- 
tistical parameter. In the Fermi limit N peaks display 
clearly. 



IV. CONCLUSIONS 

In summary, we have investigated the ground-state 
properties of ID anyon gas confined in optical lattices 
combined with a weak harmonic trap in the hard core 
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FIG. 3: Density matrix pji for trapped systems with 300 
lattice sites, Vo = 1.0 X 10 -4 i, and occupations of 50 HCAs. 
j and I range from 1 to 300. Top row: \ = 0.0 (left) and 1.0 
(right); Bottom four rows: \ = 0.2,0.4,0.6,0.8 (from top to 
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FIG. 4: (Color online) (a) Occupation of the natural orbitals; 
the modulus of the lowest natural orbital for (b) \ — 0-0, (c) 
X = 0.4, and (d) \ = 1-0. 



the one-particle Green's function we obtain the ROBDM 
and momentum distribution. It is indicated that the 
ROBDM is a complex Hermitian matrix and the mo- 
mentum distributions show properties distinct from the 
bosons and fermions. In the Bose limit and Fermi limit 
the ROBDM is real and the momentum distributions are 
symmetric about the zero momentum while anyons pop- 
ulate in the momentum space asymmetrically and would 
rather stay in some special regime with large probabil- 
ity. With the change in statistic parameter the system 
exhibits the Bose statistics, Fermi statistics and the frac- 
tional statistics in between. The anyonic system exhibits 
characteristic feature between Bose and Fermi statistics 
also in the occupation of natural orbitals and the mod- 
ulus of natural orbital shows more and more obvious os- 
cillation with the decrease in statistical parameter. 



limit using exact numerical method. With Jordan- 
Wigner transformation the hard-core anyon model is re- 
lated to polarized free spinless Fermi model and thus 
exact ground state can be constructed from that of N 
free fermions which is composed of the lowest N single- 
particle eigenstates of free fermion. Then by calculating 
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